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Abstract. We study a simplified Heisenberg spin model in order to clarify the idea of 
decoherence in closed quantum systems. For this purpose, we define a new concept: the 
decoherence function 2(i), which describes the dynamics of decoherence in the whole 
system, and which is linked with the total (von Neumann) entropy of all particles. As 
expected, decoherence is understood both as a statistical process that is caused by the 
dynamics of the system, and also as a matter of entropy. Moreover, the concept of 
decoherence time is applicable in closed systems and we have solved its behaviour in the 
Heisenberg model with respect to particle number N, density p and spatial dimension 
D in a l/r -type of potential. We have also studied the Poincare recurrences occurring 
in these types of systems: in an N — 1000 particle system the recurrence time is 
close to the order of the age of the universe. This encourages us to conclude that 
decoherence is the solution for quantum-classical problems not only in practice, but 
also in principle. 
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1. Introduction 

Decoherence is widely accepted as an explanation of how quantum correlations are 
damped out to make physical systems effectively classical. Open (infinite) quantum 
systems have been studied in great detail by many researchers [H El OH IH Ej. Their 
works are not relevant in our case, because they consider an infinite environment that 
consumes the quantum coherence irrevocably. But how about a finite environment or a 
finite system without an environment? In principle, if time is unlimited a finite system 
returns arbitrarily close to its starting position an infinite amount of times (Poincare 
recurrence). Therefore, if a finite system starts from a superposition state it begins 
to lose its coherence, but will at some moment return back to its initial superposition 
state. Within these systems, is it reasonable to talk about decoherence, and whether 
there is some kind of preferred, i.e., pointer basis as Zurek calls it [H], that is realised by 
decoherence. If the answer is "y es '\ is it somehow possible to circumvent decoherence 
in closed quantum systems? 

Our interest in closed and finite quantum systems arises from the "cosmological" 
aspects of reality. The universe has no environment [Sj and it has a finite number of 
degrees of freedom Yet decoherence is observed in our universe jH]. To understand 
decoherence, one should be able to model these critical aspects of the universe in 
decoherence studies. Previously closed quantum systems have been studied EH HI] 
using the frame of the many histories interpretation of quantum mechanics. This 
approach is, however, found to be problematic [T2*| ITH] IT1] . In this study we consider 
the off-diagonal elements of a reduced density matrix to avoid the problems of many 
histories. Another reason is that decoherence theory using reduced density matrices 
has not been studied in great detail, unlike the decoherent histories approach. A well 
established decoherence theory using reduced density matrices may clarify the concept 
of decoherence in closed systems. Moreover, these approaches are not equivalent [TU] . 

Elsewhere we have analysed in detail the conceptual problem of decoherence in 
closed and finite systems [TH]. A few major results should also be outlined here, since, 
they are the key concepts and premises. 

(i) There are two different decoherence types (similar to the different entropy types): 
the idealistic and the realistic decoherence. 

• The idealistic decoherence scheme can be applied only by those observers who 
do not interact with the universe, and who know the wave function of the 
universe and its time evolution. In a closed system, there is no idealistic 
decoherence, since the wave function of the universe remains always pure. 

• The realistic decoherence scheme is the internal view of the universe calculated 
from the wave function of the universe. It describes the events as the real 
observers that are totally correlated with the universe perceive them. To 
acquire this realistic viewpoint, one should make an effective theory from the 
wave function of the universe, e.g., to use reduced density matrices. This is 
often referred to as coarse-graining. 
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In this study, the word "decoherence" refers to the concept of realistic decoherence. 
Exceptions are mentioned. 

(ii) The possibility of recoherence does not mean that there is no decoherence. 
Decoherence is the decay of the off-diagonal elements of the reduced density matrix, 
and hence, recoherence means the growth of the off-diagonal elements of reduced 
density matrix. All finite quantum systems may experience recoherence. 

Idealistic decoherence is often referred to with the word "decoherence", since, in 
open and infinite system studies both decoherences behave similarly. The principal 
difference between idealistic and realistic decoherence can be seen only in closed systems. 

Using a reduced-density-matrix approach, we have studied a Heisenberg spin model 
in order to clarify the idea of decoherence in closed systems. We focus on one particle 
coherence (and entropy), and thereby calculate the time evolution of the system. This 
coarse-graining makes our (closed) system effectively open. Within this model, it is easy 
to sketch how decoherence is advancing in the system. The main task of this research 
is to derive functional dependences of decoherence time on relevant parameters of the 
system (particle number N, density p, potential and dimension D of the system), i.e., 
to determine how the decoherence time depends on these system variables, and why. 

This paper consists of four main sections. First, we introduce our model in section^ 
In section El we present theoretically the dynamics of a simple initial state, and define a 
decoherence function S(i) as the measure of the coherence of the whole system. Section^] 
explains the structure of our simulations, along with the main results. Section EJ is for 
discussion. 

2. The model 

We have chosen the Heisenberg spin model because it is simple enough to solve, and 
yet complicated enough to simulate properties of real quantum systems. Coupled spin 
systems are interesting from a quantum computational point of view, too. Our system, 
N interacting particles fixed in space, has no environment and, in that sense, the system 
forms a closed quantum universe. The particles are spin-| particles, and the interaction 
between them is due to their spin ^-component (analogous to the Ising model). When 
there is no coupling with the environment (i.e., no outer environment), the two spin 
states have the same energy, which is taken to be zero. Zurek jB] and Omnes ^H] have 
considered a similar, but simpler model in order to show that off-diagonal elements (i.e., 
quantum correlations) will decay in time. They labeled one particle as the system, and 
the others as an environment p, but, we study the particle system as a whole. 
The interaction Hamiltonian, 

N-l N N 

h = ky. E 9iA®< II (i) 

j=l i=j+l k=l,^j,i 
| In their model the particles that form an environment do not interact with each other. 
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describes the dynamics of the system. The interaction matrix G, where g^- = g^ gives 
the interaction strength between particles i and j. The interaction strength arises from 
the potential V; but, for formal calculations there is no need to know more about it, 
because particles are doomed to stay in one place. Fixing the positions of the particles 
is a justified assumption in decoherence studies since, in most cases, the decoherence 
time scale is the shortest time scale |2j, at least shorter than the time scale of particle 
motion. 



3. Theoretical calculations 



Let us now consider only the simplest case in order to present our method, namely 
initially a product state of superposition states 



N 



i*(o)> = n ®Ofci+fc) +h\-k)) 



k=l 



where a's and b's are normalised probability amplitudes \cik\ 2 
k = 1, . . . , N. The Schrodinger equation, 



\b k \ 



(2) 

1 for all 



ihd t \*{t)) = H\*{t)), (3) 

gives the dynamics of the system, and with the given initial condition of equation (0) 
one gets the time dependence 



\^(t)) = exp 



N-l N 

j=l i=j+l 



N 



n ® ( a k\+k) + b k \-k)) ■ 

k=l 



(4) 



The fate of the I th particle is solved by tracing over other particles, i.e., degrees of 
freedom, 

Pi = Tri,... )iV ^p, (5) 

where p = \^S(t))(^f(t)\. This is a crucial step. We make an effective theory 
of our particle system by tracing over the "uninteresting" particles (that form an 
effective environment to the particular particle of our interest), as in the mean field 
approximation. The net effect of traced-out particles is described in a simpler form and 
with lesser degrees of freedom. This results in the particular particle under consideration 
being an effectively open system. The validity of this type of coarse-graining can be 
checked by comparing the results with the entropy studies (see section fO|) . 
We thus have 



Pi = \ai\ 2 \+i)(+i\ + \k\ 2 \-i)(-i 

N 

7 * 

ai 



m n ( 

fc=l,^Z 



a fc |V i2 ^ + |6 fc | 2 e i2 ^)|+,)(-d + /i- c. 



(6) 



It is interesting that the result of equation © is the same as in reference jH], if one 
drops the index I away. In reference [0] only an interactionless environment has been 
considered, but our model counts all the interactions between particles. The result of 
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equation (JBJ) is obvious: one particle of the system will be a victim of decoherence. In 
fact, all particles separately will be victims of decoherence. Other particles act as an 
environment for the particle of interest, and the coherence of this particular particle 
is dumped (i.e., displaced) temporarily into its "environment". Now the interesting 
question concerns the decoherence of the whole particle system, as opposed to the 
decoherence of the particles separately. 

The answer can be reasoned out in the following way: let us first make our notation 
a bit lighter by denoting 

N 

z x = a x b\ J] (k|V i2 ^+|6 fc |V 2 ^). (7) 

k=l,^l 

This Zi (or its complex conjugate zf) describes the fate of the off-diagonal elements of 
I th particle. Let us put the whole system into the superposition state = bk = 
Vfc, and let the elements of the matrix G be random numbers at the interval [0, 1] §. 
The time dependencies of the off-diagonal elements for all particles of the system are 
presented in figure Note that parts of the system may come close to their starting 
level, but each recurs at different times. It is obvious that the whole system returns to 
its starting point more seldom than its parts, so the fate of all off-diagonal elements of 
the system is described by the function 

sw4ei^)I- (8) 

JV 1=1 

We have returned to a description of the whole system; but, E(t) is a quantity of an 
effective theory, and therefore our initially closed system becomes effectively open. We 
are definitely describing the system from the inside view of the system, i.e., using the 
realistic decoherence scheme. 

If S(t) achieves its initial level, the system has returned to its initial position, i.e., 
superpositions have recurred. The complement event of E(t), i.e. Y = -^J2k=i \ a kbt\ ~~ 
H(t), behaves as the statistical entropy in classical physics: it grows fast to some 
equilibrium value with certain fluctuations that depend on the number of particles in the 
system [in figureElwe have presented the behaviour of H(t)]. In open and infinite systems 
entropy and decoherence are related to each other (see, e.g., reference [2j), and therefore 
it is reasonable to assume that the same holds also in closed (and finite) systems. The 
decoherence function, H(t), is a quantity similar to the sum of all one particle (von 
Neumann) entropies of the system. The most stable (and the most probable) state is 
the one with the maximum entropy, i.e., with minimum quantum coherence. The system 
tends to reach this state and to spend most of its time in it. It should be clear that 
decoherence is a similar dynamical process as the growth of the entropy (see section [Qj) : 
therefore, recurrences in S(t) are rare and short-termed phenomena. 

§ This means that the interaction between the particles does not depend on the distances between 
particles, so, in this case, adding particles into the system means the same as increasing the density 
of the system with r-dependent interactions. It is obvious that decoherence is density-dependent, and 
thus we will consider more realistic interactions in section 0] 
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Our model shows that decoherence is advancing at different speeds at different 
parts of the system. The decoherence speed of the whole particle system differs from 
decoherence speed of the parts of the system, and it can be evaluated from the normalised 
sum of all particles, i.e., H(£). Also, the pointer basis is realised. In our model with 
a particular interaction of equation (JTJ), the pointer basis is the cr^ -basis. It is quite 
reasonable to talk about decoherence (and entropy) in a closed system. 

4. Decoherence time 

In infinite quantum systems, decoherence time is easily defined: it is the time when 
off-diagonal elements have decayed by the factor e -1 (e.g. El)- If this definition 
is straightforwardly applied in closed systems, the following problem will appear. 
Decoherence is advancing at different speeds at different parts of the system, so what is 
the decoherence time in this case? In open systems the decay of off-diagonal elements 
nicely follows the function e~' / '' rd , where is decoherence time; but, in closed systems 
the decay function fluctuates, in some cases fluctuating a lot compared to the usual 
exponential behaviour. What is the decoherence time when these fluctuations are 
present? 

We have solved this problem in the following way. First, we focus on the decoherence 
time of the whole system. Therefore, we use H(i) to solve for the decoherence time. Of 
course, it is also possible to pay attention only to some subsystem S, and solve the 
behaviour of this subsystem using 



Second, we use a least squares fit on using a function £(£) = (0.5 — c)e~ t ^ Td + c 
that allows fluctuations around the average level c. This fluctuation model is valid only 
when t is small. "Real" time averages and fluctuations are calculated after S(t) has 
stabilised. The above mentioned decay time is the decoherence time. By the way, 
the given definition yields the same results in open and infinite systems as the familiar 
form ^e~ t / Td , because in open and infinite systems c — > 0. 

4-1- Numerics and simulations 

Our main task is to try to derive the function = Td(N, p, V, D). The starting point in 
our simulations is a D-dimensional "box" whose volume is l D . N particles are placed 
randomly in this box. These particles are in fixed places, and they interact with each 
other according to the Hamiltonian (JTJ. In this paper, we concentrate on potentials of 
type 



-i m 

Zs = -J2\zi(t)\. 
m fts 



(9) 



V 



(10) 




where e = 1. Now we calculate H(t) and fit it to the function 
f (t) = (0.5 - c)e~ t/T + c. 



(11) 
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We also solve the average level of E(t) by calculating its time average over the interval 
[t 1 , t 2 ], where r d <C t 1 and Ta <C t 2 — t\. 

(Z) = F E(t)dt. (12) 



t 2 — t\ Jt-L 

This procedure is repeated U times, after which we have a statistical hunch about 
what is going on in the box. Finally, we average acquired values of t u and (E) u to 
get a statistical estimate, = jjY^i=\ T u is a function of parameters N, p, V and D. 
Standard deviations o Ti and cr(=) are also interesting. a Td is used in calculating error 
estimates of our model, and ct(h) describes fluctuations of E(t) around its average level 
(2). 

If the density under consideration is constant, then, when the number of particles 
is changed, the size of the box is changed, too: 

fN\ 1/D 

'=(7) ' (13) 

This simulation is repeated with different values of parameters (iVj, pi, V{, -D«), and 
we get results (r^, cr Td ., (H)j, C(s)J- The next task is to find the possible underlying 
functional dependence. 



4-2. Dependence on relevant parameters 

We have focused on an initial state that contains complete superpositions, so aj = 
hi — ^7= || . The number of simulations with particular parameter values is L = 100. 
Entangled particles have not (yet) been studied in relation to the case of determination 
of decoherence time, so our initial state is the state of equation ((21). The time average 
is calculated in the interval of t\ = 50 to t 2 = 100. 

In this paper, we only consider the potential (fTU|) with e = 1. The average level of 
E(t) is only a function of N, and it behaves as 

(~)(N) = Ae- BN , (14) 

where A ~ 1.27 and B ~ 0.43. The fluctuation level around (H) is the same type of 
function: a {B )(N) = Fe~ GN , where F ~ 0.03 and G ~ 0.3. 

Zurek P and Omnes [IB] have reported that theoretical fluctuations of one particle 
coherence, z(t), are ~ e~ NXa2 / 2 ~ e -o.3466A^ go ^ j g reasonaD i e that our numerical results 
for E(t) behave in a similar manner. 

The behaviour of decoherence time is more complicated, but it can also be analysed. 
The function which fits well to simulated data is 

Ti = K^ + i? ) p " (15) 

J We want to study a closed system with maximum initial coherence and minimum initial entropy and 
therefore, we set — bk — Of course, the choice of (and 6^) is arbitrary and we have tested 
that the results with afc's as random complex numbers [0, \a,k\ 2 = 1] behave similarly as in the case 
a k =b k = ^. 
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where P,Q } R } S are fitting constants which depend on dimension D, and where the 
potential scaling factor 77 is taken into account. In Table ^ we present the values of 
these fitting constants in dimensions D = 1,2,3. Dimensional dependencies may be 
extracted from these results. 

Some general results for different kinds of systems can be calculated. For D = 3 
system of iV = 100 atoms in V — 0.01 m 3 , whose interaction strength is an 
electromagnetic-type interaction (rj ~ 2.2 • 10 6 ), the decoherence time is ~ 2 • 10~ 10 s. 
We encourage experimentalists to do experiments with as good as possible isolated 
closed interacting quantum systems. 

In figure Efa2, b2, c2) differences between the dimensions are considered. The 
analysis shows clearly that H(t) describes the decoherence of the system quite well; 
especially when D > 1, it follows nicely the exponential form that is typical for open 
and infinite quantum systems. 



4-3. Recurrence 

The interesting thing to notice is that our quantum system (of N particles) is a closed 
and finite quantum system. That means, roughly speaking, that quantum correlations 
are never lost. They are only displaced, and the system may return to its initial state, if 
one waits long enough. But how much time is long enough? The answer can be reasoned 
out in the following way. Let us give an approximation about Poincare recurrence time 
Tp of our almost periodic function In reference ^Zj it has been argued that the 

period of a function of type 

F = 2J cos uj{t (16) 



is 



T Pf = . , . (17) 

mm|c<jj — LOi\ 

This is quite elementary. The same line of thinking is applicable to the products of 
cosine functions, too. So, we can give an estimate of period Tp s for the function of type 

^ N N 

S W = T?E II cos2<fct, (18) 
iV <=li=l,5« 



that is, 



(19) 



nuiiWj — ay 



where U{ = 2min|^jj — <7y/|. 

The recurrence time in our simulations can be solved, and in Table El some examples 
with various iV and D are given. The effect of p shows up only as a common factor of 
9ij in Ui. 
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Let us study the possibility for recurrences from an entropy-based point of view. 
The von Neumann entropy of a particular system that is described by the density matrix 
p is 

S(p) = -TT(p\np). (20) 

The entropy of a pure state is, of course, equal to zero. Information that lies in 
correlations between particles is given by 

N 

I = J2S(pi)-S(p), (21) 
1=1 

where index / counts subsystems, e.g., particles (see reference [IE])- In our system, 
S(p(t)) = because the system is closed and it starts from a pure state \ip(t))(ip(t)\. But 
the entropy of subsystems grows because the information, which is in the correlations 
between particles, grows. Let us now put this into a quantitative form: the total entropy 
of subsystems is 

N 

Stat = E S (Pi) = ~ E A M In Kh (22) 

1=1 i,l 

where A^ = \ ± |yl — 4(|a;| 2 |&z| 2 — |z^(t)| 2 ) are the eigenvalues of the reduced density 
matrix p\. If |a ; | 2 = |6;| 2 = |, then A^ = |± |^(t)|. In figure 121 we have plotted S to t with 
various particle numbers N. It is notable that the total entropy behaves as a mirror 
image of H(t). Moreover, the shape of the sum S(t) + * 2 S ta t{t) is similar in every 
case, even if in D = 1 there are large fluctuations present. The peak in figures |2fa3-f 3) 
indicates that, in the beginning, the entropy grows slightly faster than decoherence, but, 
after the system is relaxed, they are equal. Decoherence and entropy behave regularly; 
they are linked to each other because they have the same origin (quantum dynamics). 
Our analysis shows that E(t) is the right concept to describe decoherence in closed and 
finite systems. 

It seems clear that the system may return to its initial position, but the recurrence 
time grows fast with respect to N. From the entropy point of view, it is possible, 
in principle, for the system to get completely to the initial state again, but, it is 
thermodynamically impossible, i.e., the recurrence time is much greater than the age of 
the universe. 

5. Discussion 

We have studied this simple scenario in order to clarify the idea of decoherence in closed 
systems. We have stated that decoherence and entropy are two sides of the same coin: 
they are well defined inside (effective theory) descriptions of the closed system. This 
effective theory of a closed system is valid for observers inside the universe, i.e., they 
observe a truly closed and finite universe as effectively open. The universe is truly 
closed from an outside point of view (complete description), but the only way to access 
this view is an academic example. We leave the detailed analysis of different types of 
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observers in complete and effective theories for elsewhere. Still, decoherence is a working 
concept in closed systems also in principle, not only in practice, e.g., as Bell has argued 
[TBI EH]- Bell's argument against decoherence has been refuted (e.g., by Omnes p3J 1f- 

The possibility of recurrence makes neither concept, entropy nor decoherence, 
empty. Decoherence is a dynamical process, arising from interactions between particles, 
that diagonalises reduced density matrices in the pointer basis. This holds true in our 
model. Decoherence time can be applied and determined in closed systems as well 
as in open systems, and the function describes the fate of all off-diagonal elements 
of particles forming the system. E(t) can be linked with the entropy of the system as a 
mirror image. The particles-in-a-box example in classical statistical physics is analogous 
to decoherence: the recurrence of all off-diagonal elements in our system is a similar kind 
of phenomenon to the gas in a closed chamber going completely to the other half of the 
chamber. This is possible, but its possibility diminishes as iV grows. In the complete 
description (view from outside), the total entropy and the idealistic quantum coherence 
of a closed system are constants, but in the effective theory (view from inside) sums of 
both realistic quantities are evolving in time, as described in Sections I4.2I and I4.31 

It is quite obvious that the decoherence time behaves as ~ N~ Q and ~ p~ s , 
because decoherence is faster in bigger systems (many particles and high density) than 
in small systems. An interesting feature is that the decoherence time with respect to 
N saturates to a certain value R that is dimension-dependent. Using equation (fTKj) one 
can argue that an N ~ 43000 3-D-particle system behaves as an infinite particle system 
with an accuracy of 99%. The amount of particles needed to cause a nearly-infinite 
type of behaviour is small compared to the baryon number of our universe (~ 10 80 j7j). 
The saturation effect gives in principle a possibility to fight against decoherence: just 
make a low density system. Of course, vacuum fluctuations are present in low density 
cases, and therefore studies of the effects of second quantisation are important in fighting 
decoherence. 

With our model, it is possible study the Schrodinger cat experiment which crucially 
showed the problems in understanding quantum mechanics and the demarcation 
problem between quantum and classical [2U|. Many improvements have been made 
in quantum mechanics since the days of Schrodinger, especially concerning decoherence. 
This study covers basically the question of quantum-classical problem; we have shown 
that quantum mechanics can be applied to the whole universe, and as an effective theory 
we acquire almost a classical universe with the help of decoherence. The Schrodinger 
cat problem is more than this demarcation problem, but the presented model can be 
used in constructing a realistic enough model of a cat and its surroundings. Entangled 

% Bell's argument is mainly that if the universe starts in a pure state, it will always remain in a pure 
state, no matter how quickly the off-diagonal elements of the reduced density matrix decrease and how 
small they will become. He claims that this gives in principle possibility to make such a measurement 
that will show quantum interference. However, there is not even in principle such a measurement device 
that can perform the measurement. 

For a detailed discussion about the Bell's argument, one is encouraged to study references |15l I16| . 
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states will play a crucial role in that case. 

The "cat-in-the-box" will be in our future interests, as well as other l/r e potentials. 
Also, ideas how to circumvent decoherence are interesting. The system geometry may 
play some role in slowing down decoherence and making recurrence more probable. 
Correlations between particles should be studied as well in greater detail; here we have 
only considered an effective theory with one-particle correlations, but many-particle 
correlations may contribute to the dynamics of the system. However, our entropy studies 
seem to justify the results. There is one minor drawback in our model that is a result of 
simplifying and making an academic example, and it may affect the results: the model is, 
strictly speaking, discrete because the particles are placed (randomly, but this does not 
matter) in fixed and accurate places. Clearly, this violates the Heisenberg uncertainty 
principle and may cause considerable effects. It would be interesting to study what kind 
of results are acquired with particles as probability distributions in space. 
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Table 1. The values of fitting constants. Some constants follow (within the error 
bars which are acquired from our statistical analysis) a dimensional scaling which is 
presented on the last row. 

D P Q R S 

1 3.73 ±1.0 1.49 ±0.3 0.415 ± 0.022 1.80 ±0.05 

2 0.77 ±0.3 0.80 ±0.1 0.166 ±0.004 1.00 ±0.05 

3 0.45 ±0.1 0.55 ±0.05 0.128 ±0.003 0.67 ±0.05 
Scaling 3.5/Z) 2 1.6/D 1.80/D - 88 



Table 2. Examples of estimated recurrence times Tp- with various initial conditions. 
We have also presented standard deviations of these numerically estimated recurrence 
times. Note that Tp- is highly dependent on N and D. The unit of time is gt. 
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P 


D 


Tp s 




10 


1 


1 


5.0 


• 10 6 


1.1 


•10 8 


50 


1 


1 


4.6- 


10 10 


2.7- 


10 11 


100 


1 


1 


2.2 • 


10 12 


1.2- 


10 13 


500 


1 


1 


6.3- 


10 16 


3.6- 


10 17 


1000 


1 


1 


5.5 • 


10 17 


1.9- 


10 18 


1000 


1 


2 


2.2 • 


10 15 


4.7- 


10 15 


1000 


1 


3 


4.0- 


10 14 


8.6- 


10 14 


1000 


10 


3 


2.0- 


10 14 


4.5- 


10 14 
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Figure 1. We illustrate two cases with different N: (a), N — 4 and (b), iV = 9. For 
all particles I = 1, . . . ,N, zi(t) is plotted. We have also plotted E(t) = ^X^i \ z i(t)\ 
(a5, blO) which describes the decoherence of the whole system. Straight lines represent 
the initial states. For N = 4, recurrences are stronger and occur more often than for 
N = 9 (compare (a5) and (blO)). 
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Figure 2. These figures illustrate the relation between entropy S(t) (1) and coherence 
function (2). Various cases are considered (with p — 1): (a) N = 50, D — 1; (b) 
N = 50, D = 2; (c) N = 50, L> = 3; (d) iV = 10, = 3; (e) N = 100, 13 = 3; (f) 
iV = 500, D = 3. It seems clear that entropy (1) and coherence function (2) are mirror 
images. Also the sum 2jv 1 ln2 5'(t) + S(t) is plotted (3). 



